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Abstract 



£h ' The position of the two-proton drip line has been calculated for even-even nu- 

P' 

clei with 10 < Z < 82 in the framework of the relativistic mean-field (RMF) 

in 

theory. The current model uses the NL3 effective interaction in the mean- 
' field Lagrangian and describes pairing correlations in the Bardeen-Cooper- 

I Schrieffer (BCS) formalism. The predictions of the RMF theory are compared 

qq | with those of the Hartree-Fock+BCS approach (with effective force Skyrme 

SIII) and the finite-range droplet model (FRDM) and with the available ex- 
perimental information. 
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PACS numbers: 21.10.D, 21.10.F, 21. 60. J 



I. INTRODUCTION 

Experimental and theoretical studies of exotic nuclei with extreme isospin values are 
active areas of current research in nuclear physics. The advent of radioactive beams and the 
creation of several facilities to produce them have provided the opportunities to study the 
structure and properties of very short-lived nuclei with extreme neutron-to-proton (N/Z) 
ratios 

On the neutron-rich side, exotic phenomena include (i) the weak binding of the outermost 
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neutrons, (ii) pronounced effects of the coupling between bound states and the particle 
continuum, and (iii) regions of neutron halos with very diffuse neutron densities and major 
modifications in the shell structures. The situation is different on the proton-rich side of the 
stability valley. Here, nuclei are stabilized by the Coulomb barrier, which tends to localize 
the proton density in the nuclear interior, thereby preventing the formation of nuclei with 
large spatial extensions. 

The opportunities provided by the radioactive beam facilities make the study of the 
structure and properties of nuclei close to the proton drip line a very interesting topic from 
both experimental and theoretical points of view. Experimentally, possibilities for studying 
new decay modes such as diproton emission have opened up. Theoretical studies allow 
further tests of the various models. Of special interest is the region of sd— /p-shell proton- 
rich nuclei [|SHTT|j where two-proton ground-state radioactivity |T2|-|T5| is expected to occur. 



In particular, the region around 48 Ni is expected to contain nuclei which are two-proton 
emitters. 

In certain cases the proton drip line has been reached or even crossed experimentally. 
Systematic theoretical studies predicting the positions of the proton drip line are therefore 
important and timely []nj. In this work, the relativistic mean-field (RMF) theory is used to 
study the ground-state properties of very proton-rich, even-even nuclei with 10 < Z < 82 
and to predict the location of the two-proton drip line. 

The RMF theory ||T7|-p0f has proven to be a powerful tool to describe and predict the 
properties of nuclei. This theory provides an elegant and economical framework, in which 
properties of nuclear matter and finite nuclei, as well as the dynamics of heavy-ion collisions, 



can be calculated (for a recent review, see Ref. p0|). Compared to conventional nonrelativis 



tic approaches, relativistic models explicitly include mesonic degrees of freedom and describe 
the nucleons as Dirac particles. Moreover, the spin-orbit interaction arises naturally from 
the Dirac-Lorenz structure of the effective Lagrangian. 

In this work, the calculations are performed in the axially-deformed configuration and 
the pairing correlations are accounted in the BCS formalism. It is known that the BCS 



description of the scattering of nucleonic pairs from bound states to the positive-energy 
particle continuum produces an unphysical component in the nucleon density with the wrong 
asymptotic behavior • This effect is more pronounced for neutron-rich nuclei, for which 

the coupling to the particle continuum is particularly important. For proton- rich nuclei, 
however, the Coulomb barrier confines the protons in the interior of the nucleus. Therefore, 
the effect of the coupling to the continuum is weaker, and, for nuclei close to the proton 
drip line, the RMF+BCS approach can still be considered as a reasonable approximation 
providing sufficiently accurate solutions. Moreover, it has been shown in Ref. |Hj that the 



total energy is not affected seriously by this coupling. Of course, it is more desirable if pairing 
correlations are described in the unified framework of the Relativistic-Hartree-Bogoliubov 
(RHB) scheme [or Hartree-Fock-Bogoliubov (HFB) in the nonrelativistic approach], in which 
the nucleon densities have the correct asymptotic behavior. However, numerical codes for 
deformed RHB calculations are not yet generally available. Those appearing in published 
RHB (HFB) studies use spherical configurations [|IC|p] , |23} ^26H . On the other hand, a detailed 
study of proton-rich nuclei within the deformed HF+BCS approach with the Skyrme effective 
force SIII has been reported recently J37J . 

The current paper is the first systematic study of the proton drip-line nuclei over a 
wide range of Z values within the RMF+BCS model. In Sec. II, a brief description of the 
RMF formalism is given, while in Sec. Ill, the results of our calculations are presented and 
discussed. Ground-state properties such as binding energies, two-proton separation energies, 
proton root-mean-square (rms) radii, and deformation parameters that result from fully self- 
consistent RMF solutions have been calculated for very proton-rich nuclei near the proton 
drip line. Finally the prediction of the RMF theory for the location of the two-proton drip 
line is compared with those obtained from other theoretical models. 

Strictly speaking, the proton drip line is delineated in a Z vs. N plot by nuclei with the 
smallest positive value of the proton separation energy S\ v . To derive the global drip line, it 
is necessary to perform calculations for all nuclei, especially the odd-Z and odd-iV ones. The 
RMF calculations for these nuclei are very involved and take prohibitively long computing 



times. Therefore, this work deals only with even-even nuclei and with the two-proton drip 
line defined by nuclei with the smallest positive value of the two-proton separation energy 
5*2 P . This restriction is not too severe because it can be shown that the drip lines defined 
by Sip and S 2p are nearly parallel, except that nuclei specified by the Si p line tends, on the 
average, to have one or two fewer nucleons than those specified by the S 2p line. 



II. THE RMF FORMALISM 

In relativistic quantum hadrodynamics the nucleons, described as Dirac particles, are 
coupled to exchange mesons and photon through an effective Lagrangian. The model is based 
on the one-boson exchange description of the nucleon-nucleon interaction. The Lagrangian 
density of the model is given by pOp 

C = ip{n-d-m)4j + ^{da) 2 -U{a) 



- g a 4>a4> - g u i/ry • uip - g p ^rf ■ pTip - e^7 • A — 2 — ^ • (!) 

The Dirac spinor ip denotes the nucleon with mass m. The quantities m a , m w , and m p are 
the masses of the a meson, the uj meson, and the p meson, respectively, and g a , g^, and g p 
are the corresponding coupling constants for the mesons to the nucleon. U (a) denotes the 
nonlinear a self- interaction | |28| |, 



U{a) = l -mla 2 + ^g 2 a 3 + \g,a\ (2) 



and W v , and are field tensors |L7| 



Assuming time-reversal symmetry and charge conservation, the coupled equations of 
motion are derived from the Langrangian density (|l|). The Dirac equation for the nucleons 
is 

{-ia\7 + V(r)+{3[M + S( r )}}^= e^. (3) 
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The Klein- Gordon equations for the mesons are 

{-A + ml}a(r) = -g a p s (r) - g 2 cr 2 (r) - #f (r), 



{-A + m*}wo(r) = g w p v (r), 

(4) 

{~A + m 2 p }p (r) = g p p3(r), 
-AAo(r) = ep c (r). 

The nucleon densities act as sources, and the contributions of negative-energy states are 
neglected (no-sea approximation [jTSfl). More details on the RMF formalism can be found 
in Refs. JT/HH- 

III. NUMERICAL RESULTS AND COMMENTS 

In this work, the Dirac equation for nucleons is solved using the method of oscillator 
expansion as described in Ref. P5[ . Because most of the nuclei considered here are open- 
shell nuclei, both proton and neutron pairing correlations have been included. The BCS 
formalism was used for the pairing with constant pairing gaps obtained from the prescription 



of Ref. | 30fl . The number of oscillator shells taken into account is 12 for fermionic and 20 



for bosonic wave functions. The effective force NL3 was adopted for the calulations using a 



new version of the "axially-deformed" code |3T|. The parameter set NL3 has been derived 



recently |32| by fitting ground-state properties of ten spherical nuclei. Properties predicted 
with the NL3 effective interaction are found to be in good agreement with experimental data 



f32| , |3"3"|| for nuclei at and away from the line of (3 stability. 

The calculations have been performed for several nuclei close to the proton drip line 
for the even-even isotopic chains. In Table I the calculated total binding energies for the 
three most proton-rich isotopes close to the drip line are listed for each element with atomic 
numbers ranging from Z = 10 to Z = 82. The experimental values (in parentheses), if 



available, are also shown for comparison. With the exception of the binding energies for 
80 Zr, 100 Sn, 174 Hg, and 180 Pb, which are from Refs. p4|-|3"9|| , all other values are from the 



1995 Atomic Mass Adjustment jKj. The rms deviation between calculation and experiment 
is only 3.1 MeV. The larger differences are observed for N w Z nuclei. This observation 
might indicate that for these nuclei additional correlations should be taken into account 



4j| . In particular, proton-neutron pairing could have a strong influence on the masses. 
Proton-neutron short-range correlations are not included in our model. 
In Fig. 1, the two-proton separation energies 

S 2p (Z, N) = B(Z, N) - B{Z - 2, N) (5) 

are shown as function of the atomic number Z. In the upper panel are shown the two-proton 
separation energies for nuclei with Z = 10 — 48, while in the lower panel are shown the 
corresponding values for nuclei with Z = 48 — 82. Each curve corresponds to a given neutron 
number which changes from iV = 8 to iV = 46 (upper panel) and N = 48 to N = 94 (lower 
panel) in going from the left to the right of the figures. 

In Table II are listed (first column) the predictions of the RMF theory for the most 
proton rich even-even nuclei (with 10 < Z < 82) that are stable with respect to the two- 
proton emission, i.e., S2 P {Z,N) > 0. For comparison, the corresponding predictions of the 
HF+BCS mean-field theory (second column) with the effective force Skyrme SIII |27j and of 



the finite-range droplet model (FRDM) model |42|,|43j (third column) are also given. Finally 
in the fourth column the lightest experimentally known, proton-stable nuclei are listed for 
each even-Z element. It is seen that the predictions of the various theoretical models are in 
accordance in most of the cases. Whether such close agreement exists in the neutron-rich 
region is an open question. 

In Table III, the predictions of the RMF theory for the quadrupole deformation parameter 
/?2 are shown for all nuclei listed in the first column of Table I. It is seen that most of the nuclei 
close to the proton drip line are deformed, apart from those with magic proton (neutron) 
number, which are spherical or almost spherical. It turns out that the magic numbers 



maintain their character close to the proton drip line. In Fig. 2, the trend of the variation 
of the quadropole deformation parameter j3 2 of the most proton-rich even-even nuclei that 
are stable to two-proton emission is shown as a function of Z . 

Table III also gives the RMF predictions for the proton radii r p . Unlike the other calcu- 
lated ground-state properties, these r p values must be treated with some caution because, 
near the proton drip line, the BCS approach may not be a sufficiently good approximation 
for estimating proton radii. 

In conclusion, a systematic study of the properties of very proton-rich nuclei close to 
the drip line has been carried out. The location of the two-nucleon proton drip line has 
been predicted, which is in agreement with the predictions of other theoretical models. In 
14 of 37 cases (of even-Z elements), the proton drip line has apparently been reached in 
a variety of experiments. The existing calculations (see Table II) suggest that there are 
approximately 60 unknown isotopes of even-Z elements in the 10 < Z < 82 region that are 
proton stable. The smallness of this number reflects the increased activity in this research 
area in recent years. The number of undiscovered isotopes in the neutron-rich side is, of 
course, much larger. Calculations similar to those reported here have been carried out by 
us for over 1300 even-even nuclei on either side of the valley of stability. These results will 
be reported separately. 
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TABLES 

TABLE I. Comparison of calculated and experimental binding energies (in MeV) for some 
very proton-rich nuclei. Experimental values, where available, are displayed in parentheses. In our 
notation, 132.153 2 = 132.153 ± 0.002, 134.47 3 = 134.47 ± 0.03, etc. 



18 Ne 


134.70 


(132.153 2) 


68 Se 


572.28 


124 Nd 


1020.57 


20 Ne 


155.51 


(160.645 1) 


68 Kr 


544.35 


126 Nd 


1042.02 


22 Ne 


176.18 


(177.770 1) 


70 Kr 


575.14 


128 Sm 


1025.90 


20 Mg 


136.62 


(134.47 3) 


72 Kr 


602.92 (607.08 28) 


130 Sm 


1050.02 


22 Mg 


166.97 


(168.578 2) 


74 Sr 


605.02 


132 Sm 


1073.29 


24 Mg 


194.51 


(198.257 1) 


76 Sr 


634.86 


132 Gd 


1050.66 


22 Si 


136.94 




78 Sr 


660.08 (663.008 8) 


134 Gd 


1075.62 


24 Si 


170.61 


(172.004 20) 


78 Zr 


637.10 


136 Gd 


1098.81 


26 Si 


202.85 


(206.046 3) 


so Zr 


665.52 (669.9 15) 


136 Dy 


1075.72 


26 S 


171.17 




82 Zr 


690.59 (694.7 6) 


138 Dy 


1099.89 


28 S 


207.28 


(209.41 17) 


82 Mo 


666.70 


140 Dy 


1122.86 


30 s 


239.98 


(243.685 4) 


84 Mo 


696.05 


142 Er 


1123.66 


32 Ar 


244.56 


(246.38 5) 


86 Mo 


720.93 (725.8 5) 


144 Er 


1147.01 


34 Ar 


274.94 


(278.721 4) 


86 Ru 


698.08 


146 Er 


1171.18 


36 Ar 


302.78 


(306.716 1) 


88 Ru 


726.42 


146 Yb 


1147.13 


34 Ca 


246.29 




90 Ru 


755.03 


148 Yb 


1172.49 


36 Ca 


280.49 


(281.36 4) 


90p d 


729.27 


150 Yb 


1197.32 


38 Ca 


312.19 


(313.122 5) 


92p d 


760.26 


152 Hf 


1197.93 


40 Ti 


314.07 


(314.49 16) 


94p d 


789.17 


154 Hf 


1221.51 


42 Ti 


347.89 


(346.905 6) 


94 Cd 


762.49 


156 Hf 


1242.72 


44 Ti 


372.30 


(375.475 1) 


96 Cd 


794.21 


156-^y 


1222.58 


44 Cr 


350.43 




98 Cd 


824.87 


158 W 


1244.50 


46 Cr 


378.63 


(381.975 20) 


98 Sn 


797.11 


160 W 


1265.97 
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48 Cr 


408 


92 


(411 


462 8) 


100 Sn 


829 


46p e 


351 


34 






102 Sn 


852 


48p e 


383 


65 






i06 Xe 


874 


50 Fe 


416 


17 


(417 


70 6) 


ios Xe 


896 


50 M 


385 


20 






no Te 


918 


52 M 


418 


66 






110 Xe 


897 


54 Ni 


451 


67 


(453 


15 5) 


112 Xe 


921 


56 Zn 


452 


49 






114 Xe 


943 


58 Zn 


484 


68 


(486 


96 5) 


114 Ba 


921 


60 Zn 


510 


89 


(514 


992 11) 


116 Ba 


946 


62 Ge 


514 


11 






118 Ba 


970 




04U 


i n 

iy 


(040 




Ue 




66 Ge 


564 


71 


(569 


29 4) 


120 Ce 


974 


64 Se 


514 


40 






122 Ce 


997 


66 Se 


544 


10 






i22 Nd 


975 



94 (825.2 6) 160 Os 1244.57 

56 162 0s 1267.07 

22 164 0s 1288.71 

94 (896.70 16) 164 Pt 1267.40 

42 (919.44 6) 166 Pt 1289.23 

61 168 Pt 1310.64 

11 (921.677£) 170 Hg 1311.45 

73 172 Hg 1333.42 

37 174 Hg 1353.46 (1354.74 3) 

82 176 Pb 1354.16 

50 178 Pb 1374.40 

73 180 Pb 1394.17 (1390.65 3) 

03 

93 

49 
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TABLE II. Predictions of the RMF theory for the most proton-rich, even-even, proton-stable 
nuclei with 10 < Z < 82. Predictions of the HF+BCS mean-field theory and of the FRDM model 
are also shown. In the last column are listed the most proton-rich nuclei known experimentally. 



Calculation 
RMF+BCS 
(NL3) 



Calculation 
HF+BCS fl2 
(SIII) 



Calculation 
FRDM Kip | 



Experiment 



18 Ne 
20 Mg 

22 Si 

26 S 

32 Ar 

34 Ca 
40 Xi 

44 Cr 

46 Fe 
so Ni 

56 Zn 
62 Ge 
64 Se 

68 Kr 

74 Sr 
78 Zr 

82 Mo 
86 Ru 

90p d 

94 Cd 
98 Sn 



18 Ne 
20 Mg 
22 Si 

26 S 

32 Ar 

34 Ca 
40 Ti 

43 Cr 

46 Fe 
50 Ni 

56 Zn 
60 Ge 
64 Se 

68 Kr 

72 Sr 

76 Zr 

80 Mo 
82 Ru 

88p d 
92 Cd 

96 Sn 



18 Ne 
20 Mg 

24 Si 

28 S 

32 Ar 

36 Ca 
40 Xi 

44 Cr 

48 Fe 
50 Ni 

56 Zn 

62 Ge 

66 Se 
70 Kr 

74 Sr 
78 Zr 

84 Mo 
86 Ru 

90p d 

94 Cd 
98 Sn 



16 Ne 
20 Mg 

22 Si 

27 S 

31 Ar 

35 Ca 
39 Xi 

43 Cr 

45p e 

49 Ni 
57 Zn 
61 Ge 
66 Se 
71 Kr 
73 Sr 

79 Zr 

83 Mo 

87 Ru 
9i pd 

97 Cd 
100 Sn 
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106 Te 

110 Xe 

114 Ba 

118 Ce 
i22 Nd 

128 Sm 
132 Gd 

136 Dy 
142 Er 
146 Yb 
152 Hf 

156 W 

160 Og 
164p t 
170 Hg 
176p b 



108 Xe 

110 Xe 

114 Ba 

118 Ce 
i22 Nd 

128 Sm 
132 Gd 

136 Dy 
142 Er 

148 Yb 

152 Hf 

156 W 
162 0s 

166p t 

172 Hg 

176p b 



108 Xe 

110 Xe 

114 Ba 

118 Ce 
i22 Nd 

128 Sm 
134 Gd 

138 Dy 
144 Er 
148 Yb 
154 Hf 
158 W 

162 0s 

170p t 
174 Hg 

180p b 



106 Te 

110 Xe 
114 Ba 
121 Ce 

127 Nd 

131 Sm 
135 Gd 

l«Dy 
145 Er 

150 Yb 

154 Hf 

158 W 

162 0s 

166p t 

174 Hg 

180p b 
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TABLE III. Predictions of the RMF theory for the proton radii (r p ) and quadrupole deforma- 
tion parameters (fa) for proton-rich nuclei close to the proton drip line. 



Nucleus 




02 


Nucleus 




ft 


Nucleus 


r P 


02 


18 Ne 


2.959 


0.001 


68 Se 


4.010 


-0.285 


124 Nd 


4.854 


0.341 


20 Ne 


2.911 


0.186 


68 Kr 


4.075 


-0.274 


126 Nd 


4.862 


0.339 


22 Ne 


2.892 


0.350 


70 Kr 


4.087 


-0.310 


128 Sm 


4.905 


0.346 


20 Mg 


3.120 


0.002 


72 Kr 


4.103 


-0.358 


130 Sm 


4.911 


0.343 


22 Mg 


3.076 


0.356 


74 Sr 


4.195 


0.387 


132 Sm 


4.920 


0.341 


24 Mg 


3.021 


0.416 


76 Sr 


4.207 


0.410 


132 Gd 


4.954 


0.346 


22 Si 


3.266 


-0.001 


78 Sr 


4.213 


0.417 


134 Gd 


4.959 


0.344 


24 Si 


3.186 


0.230 


78 Zr 


4.272 


0.422 


136 Gd 


4.985 


0.359 


26 Si 


3.133 


0.320 


so Zr 


4.276 


0.437 


136 Dy 


4.998 


0.345 


26 S 


3.332 


0.001 


82 Zr 


4.205 


-0.232 


138 Dy 


5.012 


0.346 


28 S 


3.270 


0.268 


82 Mo 


4.256 


-0.230 


140 Dy 


5.017 


0.326 


30 s 


3.205 


-0.224 


84 Mo 


4.258 


-0.247 


142 Er 


5.036 


0.297 


32 Ar 


3.333 


-0.145 


86 Mo 


4.241 


0.003 


144 Er 


5.033 


0.257 


34 Ar 


3.316 


-0.176 


86 Ru 


4.308 


-0.244 


146 Er 


5.014 


-0.207 


36 Ar 


3.318 


-0.207 


88 Ru 


4.296 


0.107 


146 Yb 


5.051 


-0.251 


34 Ca 


3.393 


0.000 


90 Ru 


4.294 


0.113 


148 Yb 


5.048 


-0.207 


36 Ca 


3.375 


0.000 


90p d 


4.339 


0.109 


150 Yb 


5.049 


-0.180 


38 Ca 


3.373 


0.000 


92p d 


4.336 


0.112 


152 ffi 


5.078 


-0.163 


40 Ti 


3.524 


0.001 


94p d 


4.330 


0.071 


154 Hf 


5.062 


-0.009 


42 Ti 


3.506 


0.000 


94 Cd 


4.371 


0.071 


156 Hf 


5.089 


-0.090 


44 Ti 


3.497 


0.000 


96 Cd 


4.363 


0.003 


156 W 


5.094 


-0.006 


44 Cr 


3.607 


0.000 


98 Cd 


4.357 


0.001 


158 W 


5.117 


-0.066 


46 Cr 


3.586 


-0.004 


98 Sn 


4.394 


0.001 


160 W 


5.143 


0.110 


48 Cr 


3.603 


0.225 


100 Sn 


4.388 


0.001 


160 Os 


5.142 


0.022 



46 Fe 

48p e 
50p e 
50 M 
52 Ni 

56 Zn 
58 Zn 
60 Zn 
62 Ge 
64 Ge 
66 Ge 
64 Se 
66 Se 



3.666 
3.649 
3.655 
3.673 
3.654 
3.639 
3.810 
3.769 
3.800 
3.888 
3.904 
3.931 
3.976 
3.997 



0.003 
0.084 
0.212 
0.000 
0.001 
0.000 
0.154 

-0.001 
0.170 
0.197 
0.217 

-0.261 
0.205 

-0.265 



102 Sn 
ioe Te 

ios Te 

no Te 

110 Xe 
112 Xe 
114 Xe 
114 Ba 
116 Ba 
118 Ba 
118 Ce 
120 Ce 
122 Ce 
122 Nd 



4.411 
4.514 
4.535 
4.553 
4.600 
4.617 
4.636 
4.680 
4.717 
4.731 
4.783 
4.796 
4.805 
4.847 



0.002 
0.120 
0.142 
0.153 
0.177 
0.195 
0.221 
0.230 
0.285 
0.295 
0.315 
0.326 
0.328 
0.341 



162 0s 
164 0s 

164p t 
166p t 
168p t 
170 Hg 

172 Hg 

174 Hg 
176p b 

178p b 

180p b 



5.166 

5.189 
5.193 
5.212 
5.229 
5.254 
5.270 
5.283 
5.303 
5.313 
5.322 



-0.083 
0.106 
-0.056 
0.061 
0.066 
-0.006 
-0.001 
-0.030 
0.000 
0.001 
0.003 
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FIGURES 

FIG. 1. Calculated two-proton separation energies S^p for the N = 8 — 94 isotones as a function 
of the proton number Z. 

FIG. 2. Calculated quadrupole deformation parameters 02 of the most proton-rich, pro- 
ton-stable, even-even nuclei with proton numbers from Z = 10 to Z = 82. 
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